Polymorphic Dynamics of Microtubules 
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Starting from the hypothesis that the tubulin dimer is a conformationally bistable molecule - 
fluctuating between a curved and a straight configuration at room temperature - we develop a model 
for polymorphic dynamics of the microtubule lattice. We show that tubulin bistability consistently 
explains unusual dynamic fluctuations, the apparent length-stiffness relation of grafted microtubules 
and the curved-helical appearance of microtubules in general. Analyzing experimental data we 
conclude that taxol stabilized microtubules exist in highly cooperative yet strongly fluctuating helical 
states. When clamped by the end the microtubule undergoes an unusual zero energy motion - in its 
effect reminiscent of a limited rotational hinge. 

PACS numbers: 87.16. Ka, 82.35.Pq, 87.15.-v 



Microtubules are the stiffest cytoskeletal component 
and play versatile and indispensable roles in living cells. 
They act as cellular bones, transport roads [1| and cy- 
toplasmic stirring rods Q- Microtubules consist of el- 
ementary building blocks - the tubulin dimers - that 
polymerize head to tail into linear protofilaments (PFs). 
PFs themselves associate side by side to form the hol- 
lovif tube structure known as the microtubule (MT) . De- 
spite a long history of their biophysical study a deeper 
understanding of MT's elastic and dynamic properties 
remains elusive to this date. Besides the unusual poly- 
merization related non-equilibrium features like "tread- 
milling" and the dynamic instability there are a number 
of other experimental mysteries - in thermal equilibrium- 
that presently defy coherent explanations, most notably: 
(i) The presence of high "intrinsic curvature" [3|-[6| of 
unclear origin, also identified as a long wave-length he- 
licity [3]. (ii) In various active bending 7] or thermal 
fluctuation experiments 4] 5] 8] MTs display length de- 
pendent, even non-monotonic apparent stiffness [5|. (iii) 
They exhibit unusually slow thermal dynamics in com- 
parison with standard semiflexible filaments [5| [6| . 

The most bizarre and controversial feature (ii) has been 
the subject of much debate and some theoretical expla- 
nation attempts based on low shear stiffness modulus 
have been put forward [9|. However a careful reanal- 
ysis of clamped MT experiments. Figs 2,3 reveals two 
features not captured by these initial models: the lateral 
end-fluctuations scale as ~ L^ while the relaxation times 
scale as ~ L^. This exotic behavior naively suggests the 
presence of a limited angular hinge at the MT clamping 
point. On the other hand artifacts that could trivially 
lead to a "hinged behavior" (like loose MT attachment 
and punctual MT damage) were speciflcally excluded in 
experiments \^^. We will outline here a model based 
on internal MT dynamics explaining phenomena (i)-(iii). 
It leads us to the origin of MT helicity (i) implying (ii)- 
(iii) as most natural corollaries [10|. The two central 
assumptions of our model are as follows: (I) The tubulin 
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FIG. 1: Polymorphic Tube Model: (a) The tubulin dimer 
fluctuates between two states cr = 0, 1 (straight / curved), b) 
Tubulin switching on one MT side leads to spontaneous break- 
ing of symmetry. Combined with the built-in lattice twist MT 
forms a polymorphic helix. The polymorphic order parame- 
ter P with phase angle </> describes the distribution of tubulin 
states in the cross section, c) Polymorphic wobbling - the 
zero-energy motion of the phase angle at each cross-section. 
It is responsible for the rotation on a cone with opening angle 
a when clamped by the end. d) Defects in polymorphic order: 
Single defects (SD) have a cost proportional to their length. 
Double defects (DD) give only local energy contribution. 



dimer is a conformationally multistable entity and fluc- 
tuates between at least 2 states on experimental time 
scales. (II) There is a nearest-neighbor cooperative in- 
teraction of tubulin states along the PF axis. We are 
lead to assumptions I-II from several independent direc- 
tions: First, the experimentally observed MT helicity [3| 
implies that there is a symmetry breaking mechanism of 
individual PF's conformational properties. In analogy to 
the classic case of bacterial flagellum the existence of he- 
lices in azimuthaly symmetric bundles also necessitates a 
coo per ative longitudinal interaction along protofilaments 
[ll|[12|. Second, investigations of single protofilament 
conformations by Elie-Caille et al [l3J reveal that a single 
taxol- PF can coexist in at least 2 states with comparable 
free energy: a straight state Kpp w and a weakly curved 



state with intrinsic curvature Kpp ~ l/250nm. These 
authors also point out the apparent cooperative nature 
of straight to curved transition within single PFs. Third, 
when mechanically buckled by AFM tips tubulin dimcrs 
occasionally switch back to the initial straight conforma- 
tion 1J| . Fourth, tubulin multistability was inferred from 



the formation of stable circular MT arcs in kinesin driven 
gliding assays by Amos & Amos 115|] . Unfortunately their 
clear, seminal observations were subsequently forgotten 
for decades leading to much of the confusion about MTs 
we are witnessing today. 

Polymorphic MT Model. Starting from assumptions 
I — II we model the tubulin dimer state by a two 
state variable (T„ (s) = 0, 1 (the tubulin dimer in the 
" straight" /" curved" state, cf. Fig. la) at each lattice site 
with circumferential PF index n = 1, ...N {N = 11 — 16 
number of PFs) at longitudinal arclength centerline po- 
sition s. The total elastic -I- conformational energy can 
be written as Emt = /q (sei + etrans + einter) ds with 
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where the integration in eg; goes over the annular MT 
cross-section with ri sa 7.5nm, r2 ~ 11.5nm the inner 
and outer MT radii, with AG > the energy difference 
between the and 1 state and b w Snm the monomer 
length, J the "Ising" cooperative coupling term along 
the PF contour and with the polymorphism induced pre- 
strain Spoi oc eppCn (s) [161] where spp is the strain gen- 
erated in the curved state. The latter can be estimated 
from the switched PF curvature Kpp « (250nTO) |13j 
to be epp = dppKpp/2 « 10^^. For an isotropic Euler- 
KirchhofF beam, the actual material deformations are re- 
lated to the centerline curvature via e = — K-r with r the 
radial vector in the cross-section. 

Upon inspection it becomes clear that the phase be- 
havior (straight or curved state stability) is contained 
in the interplay of the first two terms Cei and etrans 
while the thermal dynamics is governed by the 3rd einter 
which rules over defect behavior (cf. Fig. Id). To 
understand the basic behavior we first consider a short 
MT section along which the PFs are in a uniform state 
o'n(s) — (Tn{s + h) {einter =const. cau be dropped). Fur- 
thermore we resort to the single block ansatz, i.e. at 
each cross-section there is only one continuous block of 
switched PFs of length p. This ansatz was successfully 
used by Calladine in modelling bacterial flagellin poly- 
morphic states jll| . In this approximation the energy 
density becomes 
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with the bending modulus B 
polymorphic curvature Kpoi{p) 
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Kisin(-S-p) with ki 



— j-^ — jv— . The MT phase behavior depends on the 
polymorphic-elastic competition parameter 7 = -^^^^ — 
27VAG pj^ysicaiiy^ ^ measures the ratio between poly- 
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morphic energy of tubulin switching and the elastic cost 
of the transition. For 7 < —1 all the PFs are in the 
(highly prestrained) state a = \ while for 7 > 1 all 
them are in the state cr = - both corresponding to 
a straight MT. For — 1 < 7 < 1 we have coexistence of 2 
locally (meta) stable states : straight {p — {) oi p — N) 
and curved state with p > Q. For —7 < 7 < 7 with 
7 « 0.72 the curved state is the absolute energy mini- 
mum and the straight state is only metastable. There- 
fore in this regime, the ground state of a microtubule 
bearing natural lattice twist will be helical (cf. Fig. lb). 
Assuming a stable helical state as observed in [3| we 
have p/N £ [1/4, 1/2] giving us an estimate for the ra- 
dius of curvature k~i « 9 — 14/xm. This compares favor- 
ably with an estimate of observed helices k~^ w ll/im 
from [3| . The helical stability and the magnitude of the 
protofilament curvature npp sa l/250nm [l3| with atypi- 
cal protein Young modulus Y w 1 — lOGPa, allows us also 
a simple estimate of the transition energy per monomer 
AG w -1-1.1 to +llkT. In general, the energy in Egs lHHl 
gives rise to a complex behavior and we focus on basic 
phenomena. It turns out that a most remarkable devia- 
tion from standard wormlike chain behavior comes from 
the change of polymorphic phase that we consider in the 
following. 

Polymorphic Phase Dynamics. To better under- 
stand the central phenomenon, we define at each MT 
cross-section the complex polymorphic order parameter 
P (s) = En=i e^'^^/^a,. (s) = \P (s)| e*^(«) where \P {s)\ 
denotes the "polymorphic modulus" and the "polymor- 
phic phase" (cf. Fig. lb). The polymorphic state can 
then be described by the local (complex) centerline cur- 
vature kpoi{s) = «:oe'*^*P(s) with kq = ki shnr/N and go 
the natural lattice twist that varies with PF number yjf . 
This gives rise to a helical MT shape described by the cur- 
vature \kpoi\ — Kpoi and torsion r « (j)' + qo. For large act- 
ing forces both the polymorphic phase (j) and amplitude 
\P\ will vary along the contour, however for small (ther- 
mal) perturbations the phase fluctuations will be domi- 
nant [l8| . Based on this and on the observation of stable 
helical states [3[ we will now assume \P\ = const, and 
write the total energy of the MT whose centerline deflec- 
tion is described by a complex angle 6 (s) = 9x {s)+i9y (s) 
(deflection angles in x/y direction) as follows: 

Etot = E,i {9, (/.) + Epoi (0) (5) 

The first energy term is the "wormlike-chain" bending 
contribution Eei{9,(j)) = -§ / 1^' — Kpo!| ds. The sec- 
ond term is the polymorphic phase energy Epoi (</>) = 
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FIG. 2: Effective persistence length Ip (s) as a function of 
the position from the attachment point along the MT con- 
tour. The experimental and theoretical prediction with Ib = 
25mm, A = 7.5//m, Kq^ — 18/im, qolrf, >> 1. 



-^ /n 4''^d,s with the polymorphic phase stiffness C^ — 
ksTj^ (2 + e^-'/'^'B^) which can be related to the den- 
sity of double defects with energy 2 J (cf. Fig. Id), giving 
rise to a new length scale - the polymorphic phase coher- 
ence length Z^ — C^/kBT. 

The most unusual property of a polymorphic chain is 
reflected in the rotational invariance of Epoi (0). The 
broken cylindrical to helical symmetry of the straight 
state is restored by the presence of a " Goldstone mode" 
(/) ^^ (/) -I- (/)o [lal consisting of a rotation of P by an arbi- 
trary angle 00 in the material frame (cf. Fig. Ic). This 
mode that we will call the "wobbling mode" is a fun- 
damental property of a helically polymorphic filament. 
The wobbling mode leads to dramatic effects on chain's 
fluctuations and is the clue to the resolution of myster- 
ies (i)-(iii). To see this we will first investigate the static 
properties resulting in length dependent variations of the 
persistence length. 

Persistence Length Anomalies. Among several defini- 
tions of the persistence length 
comparison with experiments 

tion persistence length l*{s) = (2/3) s'^/ {\p{s)\^) with 

p{s) = x{s) + iyis) the transverse displacement at po- 
sition s of a MT clamped at s = and (..) the sta- 
tistical average. It is easy to see from Eq. [5] that for 
small deflections, p decouples into independent elastic 
and polymorphic displacements p (s) — pei + Ppoi, such 
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FIG. 3: Longest relaxation time of a microtubule of length 
L, experimental data from [3 and theoretical prediction from 
the wobbling mode approximation and the wobbling angle 
a = 2Ko/qo ~ 7. 9° extracted from static data [30]) Fig- 2. 



the polymorphic fiuctuations induce a strong position / 
distance dependence - a behavior that could be inter- 
preted as " length dependent persistence length" . Indeed 
for /^ >> s » q^^ the persistence length displays a 
non-monotonic oscillatory behavior around a nearly lin- 
early growing average value I* (s) « |^s-f- 1^ sin {qos). 
This oscillation is related to the helical ground state while 
the linear growth l* (s) ex a^s is associated to the conical 
rotation of the clamped chain (wobbling mode), cf. Fig. 
Ic, with an angle a = 2ko'?o' • ^^^r s >> l^ the satura- 
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tion regime with a renormalized /* (00) = 1/ 1 1 

with Ipoi — 2^050 Kq is reached. The theory can now 
be compared with the experimental data [4] [5] (cf. Fig. 
2) that reveal several interesting characteristics in agree- 
ment with predictions. In particular the mean linear 
growth of /* (L) (single parameter fit /* ~ L^ gives 
5 = 1.05) and the non-monotonic Z* (L) dependence [5| 
are well captured by the theory. The linearly growing 
experimental spread of /* with L is likely linked to the 
spread of go in the MT lattice populations p/Tj. The 
large length plateau s » l^ is not reached even for 
longest MTs (~ 50pm) in agreement with coherent he- 
lices [3|. Our best comparison between theory and ex- 
periments (cf. Fig. 2) gives Ib = 25mm correspond- 
ing to y « 9GPa (proteins with Y up to 19GPa ex- 
ist [22]) and a helix wave length A « 7.5pm. This is 
close to the expected 6/xm corresponding to the twist [17| 
of the predominant 14 PF MTs fraction in the in-vitro 
MTs preparation of [Jl [5| . It turns out that Ib is larger 
than in previous studies Zb '^ 1 — &mm where however 
polymorphic fluctuations were neglected. The absence of 
the plateau also allows a lower estimate of the coherence 
length l^ > 55pm and the coupling constant J > AksT. 
Polymorphic Phase Dynamics. To describe the MT 
fluctuation dynamics we consider the total dissipa- 
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posed of an internal dissipation contribution Pint = 
h^int J ip^ds and an external hydrodynamic dissipation 
Pext = k^±J\pfds with ^^ = 47Trj/i\ni2L/r)~l/2) 
the lateral friction constant, rj the solvent viscosity, r 
and L the MT radius and length. The time evolu- 
tion equation of the phase variable (l){s,t) and elastic 
displacement pei (s,t) is given by the coupled Langevin 

-t- r^ and ^ = -^ft-- -I- r„ 



equations 4^ 



. SPdi. 



and 4^ 



with r^/p the thermal noise term. In general this dy- 
namics is highly non-linear however in the experimen- 
tally relevant regime where the behavior is dominated 
by the wobbling mode the equations simplify greatly 
and we end up with a simple diffusive behavior of the 
wobbling mode -^(j>o{t) = t^L"^/^ T^{s,t)ds with 
a friction constant given by ^tot — ^int + iext where 
^ext = 2i^Klq^'^{{l + cosLqo) - 4.sm.Lqa + qlL^ /i). For 
comparison with the experiment we compute the time 
correlation of the y deflection. A short calculation gives 
{ypoiiL,t)ypoiiL,t')) ex e-l*-*'!/"^!-^) with the relaxation 
time T {L) « L^tot/kBT. For small lengths, t (L) w 
L^int I ksT is dominated by internal dissipation while for 



large lengths r (L) 



(?o) i^. A careful analy- 



sis of the experimental data I5| reveals in fact the latter 
scaling. An independent single exponent fit gives r oc L°' 
with a = 2.9. Using the value (Ko/go) « 4. 8 x 10~^ 
from Fig. 2 and ^_l ~ 2?] with -q = IQ'^Pa ■ a [2J| we find 
the theoretical value Tth/L^ — 7.9 x lO^^s/m^ that can 
be compared with the fit of experimental data (Fig. 3) 
Tfit/L^ = 6. 25 X lO^^s/m^. The excellent agreement of 
both the exponent and the prefactor leads us again to the 
strong conclusion that in these experiments the clamped 
MT is an almost rigid helical polymorphic rotor whose 
behavior is dominated by the zero energy ("wobbling") 
mode and hydrodynamic dissipation. For very short MTs 
the linearly scaling internal dissipation dominates and 
we could measure £,int from the limit value of Tth/L, for 
L — > 0. For the available data L > 2 fim !5| this plateau- 
regime is not yet fully developed and we can only provide 
an upper estimate from the data ^i„t ^ 4 x 10~^^A^s. 

Conclusion. The MT fluctuations are well described - 
both dynamically and statically - by the bistable tubu- 
lin model and the reason for appearance of MT helices 
becomes obvious. The otherwise mysterious lateral fluc- 
tuations reflected in I* {L) ~ L^ and Tp (L) ~ L^ scaling 
are mere consequences of the "wobbling motion" of a 
polymorphic cooperatively switching helical lattice. We 
speculate that the implied conformational multistablity 
of tubulin and the allosteric interaction are not just na- 
ture's way to modulate the elastic properties of its most 
important cytoskeletal mechano-element. It could also 
be a missing piece in the puzzle of dynamic instabil- 
ity. Another intriguing possibility of using this switch 
for long range conformational signalling in vivo, could 
hardly have been overlooked by evolution. I.M.K thanks 
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Phase (p fluctuations are induced by double defects which 
carry a limited local energy cost AE = 2Jn {n being the 
number of double defects). |P| variations are induced by 
single defects having in general larger energy cost that 
grows with their end distance, cf. Fig. Id. 
This mode has a N fold symmetry however for large num- 
ber of PFs A = 11 — 15 it can be approximated as con- 
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Another more common definition, from angular correla- 
tion {cos{6 (s) — 6 {s'))} exhibiting a similarly rich be- 
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cussed elsewhere. 
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